Claims 

1 1 . A method of performing model-based optical proximity correction comprising: 

2 embedding wavefront information on a first two-dimensional complex array having a plurality 

3 of array elements and an assigned diameter; 

4 generating a phase map from said wavefront information; 

5 computing a point spread function from said phase map; and 

performing optical proximity correction calculations using said point spread function. 

1 2. The method of claim 1 wherein said step of embedding wavefront information comprises 

2 embedding simulated wavefront information from randomly generated data. 

1 3 . The method of claim 1 wherein said step of embedding wavefront information comprises 

2 embedding empirically derived wavefront data. 

1 4. The method of claim 2 comprising: 

2 forming said first two-dimensional complex array comprising a plurality of complex numbers 

3 having amplitude and phase; 

4 obtaining a complex conjugate array corresponding to said first two-dimensional complex array; 

5 arranging said first two-dimensional complex array into a symmetric complex array using said 

6 complex conjugate array; 

7 performing an analytic filtering function on said symmetric complex array; 

8 applying a power law function to said symmetric complex array; 

9 computing said wavefront by performing a Fourier Transform on said symmetric complex array, 

I o resulting in a transformed array; 

I I obtaining a circular core of said transformed array of a diameter equal to said first array's 

12 assigned diameter, and converting to zero coordinates outside said diameter; and 

1 3 normalizing amplitudes of said transformed array. 

1 5. The method of claim 4 wherein said step of forming said first two-dimensional complex 
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2 array further comprises putting said array elements in phasor form with a phase within a range of 

3 2n. 

1 6. The method of claim 5 including said phasor having an amplitude in the range of 0 to a 

2 predetermined small fraction less than 1 . 

1 7 . The method of claim 5 wherein said first two-dimensional complex array size is represented 

2 by the log in base two of said first two-dimensional complex array's assigned diameter. 

1 8. The method of claim 4 wherein said complex numbers are centered within said first two- 

2 dimensional complex array. 

1 9. The method of claim 4 wherein said step of arranging said first two-dimensional complex 

2 array into said symmetric complex array further includes equating each amplitude of said first two- 

3 dimensional complex array with a corresponding amplitude of said complex conjugate array. 

1 1 0. The method of claim 4 wherein said step of arranging said first two-dimensional complex 

2 array into said symmetric complex array further includes equating each phase of said first two- 

3 dimensional complex array coordinate pair with the inverse of a corresponding phase of said 

4 complex conjugate array. 

1 11. The method of claim 4 wherein said step of performing an analytic filtering function 

2 includes multiplying said amplitude of each of said complex number by an apodizing function. 

1 12. The method of claim 1 1 further including filtering out low spatial frequency aberrations. 

1 13. The method of claim 4 wherein said step of applying a power law function on said 

2 symmetric complex array comprises multiplying said amplitude of each of said complex number 

3 by the sum of the squares of each of said coordinate pair, said sum raised to a power. 
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1 14. The method of claim 13 wherein said power comprises a supplied user function in the form 

2 of y/2 with y in the range of 1 to 3. 

1 1 5. The method of claim 4 wherein said Fourier Transform is a two-dimensional Fast Fourier 

2 Transform. 

1 16. The method of claim 15 including assigning each phase of said complex number to zero 

2 after said transform. 

1 17. The method of claim 4 wherein said step of normalizing said transformed array's amplitudes 

2 comprises: 

3 subtracting an average value of said wavefront from each array element; and 

4 replacing each of said elements by an analytical expression, which is a function of both a square 

5 root of intrinsic flare and a radius of said wavefront. 

1 18. The method of claim 4 including: shuffling said elements such that said element in said first 

2 array's center is shifted to said first array's bottom-left corner. 

1 19. The method of claim 1 8 further including: 

2 ignoring a first row and first column of said first array; 

3 dividing said first array's remaining portion into four quadrants; 

4 exchanging a first quadrant with a second quadrant, and exchanging a third quadrant with a 

5 fourth quadrant; 

6 exchanging said first and said fourth quadrants; 

7 exchanging said second and said third quadrants; 

8 converting said elements into real and imaginary form; 

9 unshuffling said elements such that said first array is again centered following said Fourier 
10 Transform; and 
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1 1 converting said first array into phasor form. 

1 20. The method of claim 3 further comprising: 

2 reading said empirically derived wavefront data in a row major order by substituting zero for 

3 ignored data; and 

4 centering said empirically derived wavefront data. 

1 21. The method of claim 20 further comprising: 

2 obtaining a circle enclosing said empirically derived wavefront data; and 

3 embedding and centering said circle within a square array such that said circle has a diameter 

4 represented by an equal number of rows and columns of said array expressed as a value 

5 equal to 2 raised to the power of the log of the sum of said diameter plus one. 

1 22. The method of claim 1 9 further comprising: 

2 exponentiating empirically derived wavefront information resulting in an array of complex 

3 numbers in phasor form; 

4 embedding said exponentiated wavefront in a larger array having a center such that if a linear 

5 dimension of said wavefront array is given by log diameter p, then a corresponding linear 

6 dimension of said larger array is given by log diameter q, where q is at least p+3; 

7 shuffling said larger array to move said center to a bottom-left quadrant; 

8 performing an inverse Fast Fourier Transform on said larger array to obtain a transformed array; 

9 unshuffling said transformed array to move said complex numbers from said bottom-left 

1 0 quadrant back to said center; 

1 1 normalizing the magnitude of the unshuffled array; 

1 2 choosing amplitude values of said normalized array and creating a phase map; 

1 3 using said phase map to compute a point spread function array; and 

14 scaling said point spread function array from pixel dimension to real dimension. 

1 23 . The method of claim 22 wherein said array of complex numbers further comprises phase 
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2 values which are the same as corresponding amplitude values. 

1 24. The method of claim 22 including amplitude values assigned to unity. 

2 25. The method of claim 22 further wherein said large array includes a number of rows and 

3 columns calculated by a value equal to 1 added to 2 raised to the p power. 

1 26. A method of performing model based optical proximity correction on a lithographic mask 

2 pattern incorporating phase maps comprising: 

3 incorporating a point spread function array in real dimension within a set of convolution 

4 kernels; and 

5 computing an aerial image with aberrations using said set of convolution kernels. 

1 27. A program storage device readable by a machine, tangibly embodying a program of 

2 instructions executable by the machine to perform method steps for performing model-based optical 

3 proximity correction, said method steps comprising: 

4 embedding wavefront information on a first two-dimensional complex array having a plurality 

5 of array elements and an assigned diameter; 

6 generating a phase map from said wavefront information; 

7 computing a point spread function from said phase map; and 

8 performing optical proximity correction calculations using said point spread function. 

1 28. The program storage device of claim 27 wherein said step of embedding wavefront 

2 information comprises embedding simulated wavefront information from randomly generated data. 

1 29. The program storage device of claim 27 wherein said step of embedding wavefront 

2 information comprises embedding empirically derived wavefront data. 

1 30. The program storage device of claim 28 comprising: 
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2 forming said first two-dimensional complex array comprising a plurality of complex numbers 

3 having amplitude and phase; 

4 obtaining a complex conjugate array corresponding to said first two-dimensional complex array; 

5 arranging said first two-dimensional complex array into a symmetric complex array using said 

6 complex conjugate array; 

7 performing an analytic filtering function on said symmetric complex array; 

8 applying a power law function to said symmetric complex array; 

9 computing said wavefront by performing a Fourier Transform on said symmetric complex array, 

1 0 resulting in a transformed array; 

1 1 obtaining a circular core of said transformed array of a diameter equal to said first array's 

12 assigned diameter, and converting to zero coordinates outside said diameter; and 

1 3 normalizing amplitudes of said transformed array. 
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